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We consider the existence in arbitrary finite dimensions d of a POVM comprised of d 2 rank-one 
operators all of whose operator inner products are equal. Such a set is called a "symmetric, in- 
formationally complete" POVM (SIC-POVM) and is equivalent to a set of d 2 equiangular lines in 
C d . SIC-POVMs are relevant for quantum state tomography, quantum cryptography, and foun- 
dational issues in quantum mechanics. We construct SIC-POVMs in dimensions two, three, and 
four. We further conjecture that a particular kind of group-covariant SIC-POVM exists in arbitrary 
dimensions, providing numerical results up to dimension 45 to bolster this claim. 



I. INTRODUCTION 



In quantum theory, measurements are represented by positive operator valued measures (POVMs). A POVM is 
termed informationally complete if its statistics determine completely the quantum state on which the measurement 
is carried out [1-3, 16]. In order to be maximally efficient at determining the state, such a measurement should 
also be rank-one; i.e., the measurement operators or POVM elements should be positive multiples of projectors onto 
pure states, in which case each POVM element corresponds uniquely (up to a phase) to a subnormalized vector in 
C d . A particularly appealing and potentially useful measurement is one which is symmetric, meaning all pairwise 
inner products between the POVM elements are equal. Such a POVM is a "symmetric, informationally complete 
positive operator-valued measure," or SIC-POVM for short. The set of vectors comprising a SIC-POVM has also 
been studied in a very different context, where it has a different name: it is a set of d 2 equiangular lines in C d , 
first studied by Lemmens and Seidel [4] and subsequently by many others [5-12]. In quantum information theory 
such measurements are relevant to quantum state tomography [13], quantum cryptography [14], and to foundational 
studies [15] where they would make for a particularly interesting "standard quantum measurement" . The outstanding 
question we address in this paper is whether SIC-POVMs exist in any finite dimension. 

We conjecture that SIC-POVMs exist in all finite dimensions and, moreover, that there exists in all hnite dimensions 
a SIC-POVM that is covariant under a standard representation of x Z^. To state the conjecture, let us formalize 
the definition of a SIC-POVM. The simplest definition is that a SIC-POVM P is a set of d 2 normalized vectors \<t>k) 
in C d satisfying 

More precisely, the POVM elements of P are the subnormalized projectors \4>k)(4>k\/d = Hk/d, which have pairwise 
Hilbert-Schmidt inner product (IP, II fe )/d 2 = Tr[Il]lT fc ]/rf 2 = l/d 2 (d+ 1) for j ^ k. It turns out that the other 
properties of a SIC-POVM, i.e., completeness and informational completeness, follow from Eq. (1), as we show in 
section 2. 

We can now state our conjecture. 

Conjecture 1 For any dimension d <E N, let {|fe)}t.Zj be an orthonormal basis for C d , and define 

d-i 

uj = exp{2ni/d) , D jk = uj lk/2 ^ (J m \k © m){m\ , (2) 

m=0 

where © denotes addition modulo d. Then there exists a normalized \<j>) £ C d such that the set {Djk\(f)}j fe=1 is a 
SIC-POVM P. 

Analytic solutions are known for g? = 2,3,8 [8], and to this list we add d — 4. Additionally, computer calculations 
reveal numerical solutions (with an accuracy better than 1 part in 10 s ) in dimensions up to 45, some derived using 
the aforementioned group, but others using other suitable groups. These results are detailed herein according to the 
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following plan. Section 2 states the problem in the language of frame theory and derives a connection to the problem 
of finding spherical t-designs. Section 3 specializes to the group-covariant case and explains why our conjecture might 
be generally true. Section 4 presents our analytic solutions for d=2, 3,4, and section 5 the method of obtaining the 
numerical results. Finally in section 6 we discuss possible approaches to a general proof, as well as related open 
questions. 

II. FRAMES AND SPHERICAL t-DESIGNS 

The concepts of frame theory provide a simple and elegant means of putting our problem in a general setting, 
for a SIC-POVM is a particular kind of frame. Frames are a generalization of basis sets, with the requirements of 
orthogonality and normalization relaxed. For a finite-dimensional vector space TC, a collection of vectors \ipk) <= W is 
a frame if there exist constants < a < b < oo such that 

a\m\ 2 <Y.\^k)\ 2 <b\m\ 2 (3) 

k 

for all |£) e H. Any collection of vectors is a frame in the subspace spanned by the vectors. The constants a and b 
are called the frame bounds, and if a = b, the frame is said to be tight. The frame operator is the positive operator 

S = £>k>Wfcl- (4) 
k 

It should be immediately clear that for a tight frame S = al. This tight-frame condition is equivalent to the com- 
pleteness condition for the corresponding POVM elements \ipk}(ipk\/a, and thus rank-one POVMs and tight frames 
are the same mathematical object. 

Now let § d C C d be the subset consisting of vectors that have unit norm. Any frame can be rewritten in terms 
of the corresponding normalized vectors, but tightness is not preserved under this transformation. For a frame 
{\tpk) € made up of normalized vectors, the quantity 

Tr[S 2 ]=]T|<^# fc >| 2 (5) 

j,k 

is called the frame potential. We consider only frames made up of normalized vectors throughout the following. 
A useful theorem due to Benedetto and Fickus [18] states the following. 

Theorem 1 (Benedetto-Fickus) Given any d and n, let {IV'fe) G § d }/J =1 be a set of normalized vectors with frame 
operator S. Then 

Tr[S 2 ] > max(n,n 2 /d) . (6) 

Furthermore, the bound is achieved if and only if {\"4>k)} consists of orthonormal vectors, when n < d, or is a tight 
frame, when n > d. 

Proof The proof is so simple that we can include it here for completeness. Denoting the ordered eigenvalues of S 
by Ai > A 2 > . . . > Arf, we first note that the number of nonzero eigenvalues is at most q = min(n, d). Thus we have 

q q 
Tr[S] = n = J2 X k and Tr^ 2 ] = \\ . (7) 

fc=l k=l 

Minimizing Tr[S" 2 ] subject to the constraint Tr[S] = n gives the inequality. Equality holds if and only if A^ = n/q, 
k = l,...,q. Thus for n < d, S is a projector onto an n-dimensional subspace, implying that the vectors \ipk) are 
orthogonal, and if n > d, S = (n/d)I, implying that the set {IV'fc)} is a tight frame. □ 

Since the frame potential of P is TrfS* 2 ] = d 3 =n 2 /d, the theorem establishes immediately that P is a tight frame, 
whence P is also a POVM. 

For P to be informationally complete, the d 2 operators 11^ = \<ftk)(<ftk\ must be linearly independent so that 
they span the space of operators. The linear independence follows from considering the rank of their Gram matrix 
(Ilj,rifc) = Tr[njn fe ] = (dSj k + l)/(d + l), which being circulant (each row is a cyclic shift of the previous row), 
has eigenvalues given by the Fourier transform of one of the rows. A simple calculation reveals that due to the 
combination of constant term and Kronecker delta, the eigenvalues are exactly the same as the values in any row. 
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Since no eigenvalues are zero, the Gram matrix has full rank, the projection operators Ilk are linearly independent, 
and P is informationally complete. 

Since every rank-one POVM is a tight frame, a SIC-POVM P is clearly something more. To fully elucidate the 
properties and applications of P, we need to introduce spherical t-designs. Building on the result of Benedetto and 
Fickus, we can establish a connection between frames and spherical t-designs applicable to the SIC-POVM problem. 

A spherical t-design is a set of n normalized vectors {\<frk) € S d } such that the average value of any t-th order 
polynomial ft(ip) over the set {|</>fc)} is equal to the average of /t(V') over all normalized vectors Note that if a 
set is a t-design, it is also an s-design for all s < t, since an s-th order polynomial is also a t-th order polynomial. 
Spherical t-designs were originally developed as subsets of the real sphere S d ; here we apply the concept to the set S d . 

Let H = C d , H t be the t-fold tensor product of such spaces, and S t be the symmetric subspace of H t , and consider 
a function f t :H^C defined as 



m = 



e H , 



(8) 



where the choice of ft is equivalent to a choice of a symmetric operator F t G B{St). Such a function is a t-th order 
polynomial function on H. We can decompose F t into a sum of product operators, i.e., Ft=J2k ®j=i Aj';fci tnus an y 
such function can be decomposed into monomial terms like 



**)=n< 

3=1 



Without loss of generality, we can restrict our attention to such monomial functions and rewrite them as 



/ t (vo=n Tr [^#)<v>i 



3 = 1 



= Tr 




\3 = 1 



(9) 



(10) 



Since the set {|0fc)} is a t-design if and only if the average of any f t over {| is equal to its average over all 
!>) G § d , we are led to compute the average of an arbitrary monomial term: 



(ft) 



d^Tr 




\3 = 1 



= Tr 







= Tr 









(11) 



Hence we focus on finding K t , since it effectively takes the average of f t . A spherical t-design is then a set of vectors 
for which 



St = j2\*k)(*k\ = nK t' i*k> 



(12) 



fe=i 



Note that S t is the t-fold tensor-product analog of the frame operator S. 

To find the operator K t , note that K t has support only on the symmetric subspace S t . Further, because K t is 
invariant under any U m for U € SU(d), we conclude that K t oc n sym , the projector onto S t . (Recall that <S t is an 
irreducible invariant subspace of the group consisting of the operators U® 1 .) Finally, to determine the constant of 
proportionality, we consider the average of the trivial function f t (ip)~l. Equation (11) then becomes Tr[.Kt] = l, and 
since S t has dimension ('"jjl^ 1 ), we have 



_ t!(d-l)l 

t_ (t + d-iy. sym 



(13) 



For t = 1, we see that a 1-design is a tight frame made up of normalized vectors and, hence, also a POVM made up 
of equally weighted rank-one projectors. 

Equation (12) now says that the set {|<?ife)} is a t-design if and only if the set {|$|.)} is a tight frame on St, whence 
we can apply Theorem 1 to obtain the following result. 

Theorem 2 A set of normalized vectors {\<f>k) <= w ^ n > (^l^ 1 ) forms a spherical t-design if and only if 



nH\ {d-iy. 
(t+d-iy. 



(14) 



Furthermore, this value is the global minimum ofTr[St]- 
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This theorem links the spherical i-design property with the minimization of the t-th frame potential, Tr^S 2 ]. 
Immediately we can infer that every SIC-POVM is a 2-design since Tr[Sf] = £\ fe \((f>j\(f>k)\ 4 = 2d 3 /(d+ 1), the 
required value for a 2-design. The converse is also true, namely, every 2-design with n = d 2 elements is a SIC-POVM. 
To show this, let Xjk = |(<^j|<^fc)| 2 , j ^ k, and interpret these Xjk as coordinates in R d ( d Using the values of the 
frame potentials for a 2-design with n = d 2 elements, we can write 

^A jfc =Tr(^)-^ = ^^ and £ X% = Tr(S 2 ) - d 2 = . (15) 

The first equation describes a plane and the second a sphere. They intersect at the single tangent point Xjk = 1/ (d+l), 
thus showing that all 2-designs with d 2 elements are SIC-POVMs. These considerations make clear that the crucial 
distinguishing property of a SIC-POVM is that it is also a 2-design. Moreover, this ensures that minimizing the 
second frame potential, as we do in the numerical work reported in section 4, yields vectors that do indeed form a 
SIC-POVM. 

Furthermore, d 2 is the smallest number of elements a 2-design can have, so a SIC-POVM is a minimal 2-design 
(see also [8]). Theorem 2 does not provide the minimum number of states, n m ; n , for a i-design in d dimensions, but 
we can establish lower bounds. For the t = 2 case, consider again the steps leading to the definition of the operator 
K t . In carrying out the average of the function / 2 , we could have written 

(/ 2 ) = Tr[A 2 j d^ IVXVI^IVXVl] - r ^[A 2 g{A 1 )] (16) 

and thus considered the superoperator Q : B(C d ) -» B(C d ). Here G{UAW) = UG(A)W for any U e SU(d), so 
by Schur's lemma, Q is some linear combination of projectors onto the invariant subspaces of U acting on B(C d ). 
These invariant subspaces are (i) the (d 2 — l)-dimcnsional subspace of traceless operators and (ii) the one-dimensional 
subspace spanned by the identity operator /. Thus we can write Q = al + bl where T{A) = A and 1(A) = Tr[A]I (i.e., 
1 is the identity superoperator, and I projects onto the identity operator). To find a and 6, we first let A\ =A 2 =I, 
which gives the function / 2 (V>) = 1, so that Eq. (16) yields d(a + bd) — 1. Next we consider A\ = A^ = \4>)((f)\, for 
which (/ 2 ) = fdip |(0|V')| 4 = a + b. We can use Eqs. (11) and (13) to show that (/ 2 ) =2/d(d + 1); combined with the 
previous result, this implies a = b= l/d(d+ 1). Therefore, Q has no null subspace, must be rank-<i 2 , and cannot be 
constructed from less than d 2 linearly independent rank-one superoperators. Similar arguments can be applied to all 
spherical ^-designs. By similar rearrangements, we can make several different types of operators K' t , and the rank of 
each serves as a lower bound on the number of vectors required to comprise a i-design. 



III. GROUP COVARIANCE 



Our results are obtained by considering group-covariant sets, so it is appropriate to specialize to this case. The 
SIC-POVM P is group covariant if there exists a group G with a d-dimcnsional projective unitary representation 
{U g } such that (i) P is invariant under any U g , i.e., for any \cf>j) e P and any U g , U g \(j)j) G P (up to a phase), and 
(ii) {U g } acts transitively on P, i.e., for any \<j>j), \4>k) € P, there exists U g such that U g \<j)j) = \4>k) (also up to a 
phase). Assuming group covariance simplifies the search for SIC-POVMs. We simply search for a fiducial vector such 
that P — {U g \<p)} is a SIC-POVM (note that the transitivity property implies that the order of G must be at least 
d 2 ). To do this, we use groups such that {U g \(f))} is a 1-design, i.e., a POVM, for any normalized vector \(f>), and then 
we search for a particular vector \<j>) such that |(</>|£/g|</>)| 2 = l/(<i+ 1) for all g ^ e. All other inner products are then 
guaranteed to have this value due to the group action. 

We suspect the case of group covariance to be general for the following reason. Consider the map a : S d — > B(C d ) 
that takes a normalized vector to the corresponding projector, i.e., a(\<pj)) = \<t>j){^>j\- Now consider the operators 



d-i^w-i) ■ (17) 

Being both traceless and Hermitian, these operators lie in a subspace of B(C d ) that is isomorphic to M. d _1 ; indeed, 
since (<jj,<Tj) — 1, they all lie on the unit sphere in R d _1 . This sphere is a generalization of the Bloch sphere for 
two-dimensional systems, the difference being that for d > 2, not all operators on the sphere are images of vectors in 
§ d under the map a. From the SIC-POVM condition (1), one finds immediately that 

(^•.^) = -^TI V -^ fc - ( 18 ) 
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This is the condition for the d 2 operators {<Jj} to form a regular simplex in M. d _1 , whose automorphism group is the 
permutation group S^. Given this result, some group covariance seems natural. One is tempted to think that from 
here it is a simple matter to establish the existence of the set P. This is not the case, however, as working in the 
operator space obscures the very difficult task of determining when a given operator is the image of some element of 
§ d under the map a. In the same vein, most of the elements of the permutation group cannot be represented in this 
framework as unitary transformations of C d ; thus, while we know that any G satisfying the conditions above must be 
a subgroup of S,p, it is not obvious which subgroups are candidates. 

The outstanding choice for G is the group Zd x Z^, as described in the conjecture. We note in passing that using 
this group to find a SIC-POVM makes P a Gabor or Weyl-Hcisenbcrg frame, this being the definition of such a frame 
[24]; such a SIC-POVM P is useful in defining finite-dimensional analogs of the familiar P and Q quasidistributions 
of infinite dimensions. The group's usefulness here stems from the fact that, for any normalized G § d , 

S^ = Y,D ik \i>){i>\D} k =dI , (19) 

jk 

a fact readily checked by direct calculation. 

The property in Eq. (19) of producing a 1-design for any input state is quite general. The following argument is 
adapted from Proposition 3 of [17]. Any set of d 2 orthogonal unitary operators Tj, thus satisfying Tr[TjT k ] = dSjk, 
is a complete set for expanding operators in B(C d ); the unitary operators {Dj k } are one example of operators that 
satisfy this orthogonality condition. It is a simple matter to turn the completeness relation into Y^ fc T k CT k = dTr[C]I 
for any operator C. Simply consider the inner product of two arbitrary operators A and B. The completeness relation 
means that 

(A,B) = ^(A,T k )(T k ,B). (20) 
a k 

Setting A= \(j>\){4>2\ and B= l^iX^I (which we can do without loss of generality because such outer products span 
B(C d )), we find 

(faWiXHfa) = \^(<h\Tk\<h)(HTi\ih) , (2i) 

k 

from which it follows that 7fc|(/> 2 )(V>2|7fc = d{ip2\4>2), whence the result follows. 

Thus the property of producing a 1-design regardless of the fiducial \<j> ) is common to all groups of size d 2 whose 
representation operators are a complete, orthogonal set. Such groups were introduced by Knill in connection with 
quantum error-correcting codes and are called "nice error bases" or unitary error bases [22, 23]. Klappenecker and 
Rotteler have kindly detailed all such nice error bases up to dimension 10, so we can apply them to the problem at 
hand [19-21]. Only the nice error bases associated with the group Zd x Zd exist in every dimension, thus accounting 
for our focus on this group. 



IV. ANALYTIC SIC-POVMS 



Here we concentrate specifically on using the group Zd x Zd- Fixing the representation operators Dj k of this 
group we can determine the set of fiducial vectors that under the group action make a SIC-POVM. From this we can 
determine also the number of distinct SIC-POVMs generated by our fixed representation. In three dimensions there 
are an uncountably infinite number of such covariant SIC-POVMs, but in two dimensions there are just two, and in 
four, 16. We write the fiducial state as 

\d>) =J2r k e ie "\k) , (22) 
k 

where we can, of course, immediately choose O = 0. 



A. d= 2 

The two solutions, represented as column vectors in the standard basis, are 

fj_/ V3 + V3 \ J_/ -V3- V3 \1 
{V6\ e-/ 4 73 / ' V6 \ e-/ 4 y3+73 J J 
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These have a simple interpretation on the Bloch sphere, where the nontrivial group operators are simply rotations by 
7r about the x, y, and z axes, respectively. Then the Bloch vectors of the two fiducial states are ±(1, 1, l)/\/3, and 
the two SIC-POVM states thus formed are regular tetrahedra, each one related to the other by inversion of the Bloch 
vectors. 



For r satisfying 1/V2 < r < y/2/3, define 



B. d = 3 



1 1 



r±(ro) = |ro±| A / 2 - 3r o • 



(24) 



Hence < r_ < 1/V6 < r + < 1/V2 < r < \/2/3. The complete set of fiducial states, represented as column vectors 
in the standard basis, is then 



js 1 



( plus all vectors formed 



»e 2 / ' \ by permuting of elements 




plus all vectors formed 
by permuting of elements 



< 1 < 2tt } . 



(25) 



C. d = 4 



Now let 



along with 



1 - l/VE , ,- 

ro= r —L- , r 1 = (V2-l) 
2V2-V2 



r± = iyi + l/V5±^l/5 + l/V5 



2 k ■ 2 

a — arccos — -j^=^= , b = arcsin —= , 

V 5 + V5 V5 

and define the set 

n = \ ((-l) m (a/2 + 6/4) + 7r(m + 2n + 7j + l)/4, 7r(2fc + l)/2, 



(-l) m (-a/2 + 6/4) + 7r(m + 2n + 3j + 4fc + l)/4 



j, fc, m = 0, 1 and n = 0, . . . , 3 



The complete set of fiducial states, represented as column vectors in the standard basis, is now 

/ plus all vectors formed i (000) Q 
' I by cycling of elements ' ^ + ' 1 ' ~ ' 





( ro 


\ 


( ro 


\ 




r+e t9 + 




< 


ri e idl 










\ r-e %e - 


) 


\ r + e ie + 


J 



(26) 



(27) 



(28) 



(29) 



V. NUMERICAL SIC-POVMS 

Because analytic solutions to the SIC-POVM condition (1) are so few, our conjectures are based almost entirely 
on numerical evidence (even the d = 4 solution was originally inspired by close examination of numerical solutions) . 
To find numerical instances of P, we simply minimize the second frame potential Tr^f] over sets of d 2 normalized 
vectors generated by a representation of x "L d from a vector \<j>). It is also possible to vary independently the d 2 
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elements of P, but this is much less efficient; taking advantage of the group-covariance conjecture permits us to search 
a space of 0(d) complex parameters instead of 0(d 3 ) complex parameters. 

The quantity that we minimize, J2j k \{ < t > Wjk\ < t>}\' 1 1 is proportional to the frame potential because of the group 
covariance. Since it is a quartic function of \<p), we have to use numerical methods to minimize it, using cither 
Mathcmatica (simpler) or C++ (much faster). The method used is an adaptive conjugate gradient method; this has 
the advantage of converging with exponential rapidity to a local minimum, but the disadvantage of being insensitive 
to global conditions. As a result, the most time- intensive portion of the computation by far is identifying one of the 
global minima among the many local minima. 

Once the correct minimum is located, we quickly obtain P such that Eq. (1) is satisfied to an accuracy of 10~ 8 . The 
sole exception to this rule is d = 3 (where an exact analytic solution is known): in d = 3 there exists a continuously 
infinite family of solutions, and this degeneracy makes numerical solution difficult. For every dimension between d = 5 
and d = 45, however, we have found Z^ x Z^-covariant solutions to within machine precision [25]. 

Additionally, in small dimensions, one can attempt an exhaustive search for all possible Z<j x Z^-covariant SIC- 
POVMs, by simply running the minimization many times with differing presumptive fiducial states, tabulating all the 
while the distinct SIC-POVM fiducial states found. Table 1 lists the results for the number of distinct SIC-POVMs. 



d 


#(SIC-POVMs) 


2 


2 


3 


oo 


4 


16 


5 


80 


6 


96 


7 


336 



TABLE I: Number of SIC-POVM sets generated by a fixed representation of the group Z d x 7h d in dimensions two through 
seven. The infinity in dimension three is uncountable. 

Finally, we have tested some of the other nice error bases tabulated by Klappenecker and Rotteler. These are 
also easy handled, and although not all groups were tested, at least four groups were found to generate SIC-POVM 
sets. In the notation of the library of small groups used by GAP3, GAP4, and MAGMA, these groups are G(36,ll), 
G(36,14), G(64,8), and G(81,9). Each of these solutions has an accuracy of 10~ 15 in the individual vector inner 
products. Perhaps surprisingly, many of the tabulated groups do not seem to yield group-covariant SIC-POVMs. 

VI. ODDS AND ENDS 

A rigorous proof of existence of SIC-POVMs in all finite dimensions seems tantalizingly close, yet remains somehow 
distant. Although the numerical evidence makes very clear the relevance of the group x Z,;, this is not definitively 
established. Given the apparent importance of Z^ x "L d , it would seem to be just a short step to some general form 
for an operator whose eigenvectors could be a fiducial state, but a proof by this method has not been forthcoming. 
For instance, in three dimensions the Fourier transform operator has an eigenvector that is a fiducial state (the one 
associated with the eigenvalue i), but this does not hold in general. In five dimensions a fiducial vector can be found 
among the degenerate eigenvectors of a particular Z3 subgroup of the normalizer of Z^ x Z^ in SU(d), but there 
is no such subgroup at all in the normalizer for dimension seven. The group-theoretic structure of SIC-POVMs is 
exceedingly rich, however, and ongoing efforts to understand the full automorphism group of a SIC-POVM might 
yield insights into operators that yield fiducial states. We hope that by establishing the framework and providing 
motivating numerical results in this paper, a proof might be completed. 

Finally, regardless of whether their existence can be proved rigorously, SIC-POVMs appear to exist in many di- 
mensions. They might be of use in many areas of quantum information theory — signal ensembles, quasidistributions 
on discrete phase space, error-correcting codes, and quantum measurement, just to name a few. In addition, they 
appear to be connected to a number of interesting mathematical problems, including spherical codes and the putative 
existence of mutually unbiased bases. 
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